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We study the possibility of constraining the energy dependence of cosmological birefrin- 
gence by using CMB polarization data. We consider four possible behaviors, characteristic of 
different theoretical scenarios: energy-independent birefringence motivated by Chern-Simons 
interactions of the electromagnetic field, linear energy dependence motivated by a 'Weyl' in- 
teraction of the electromagnetic field, quadratic energy dependence, motivated by quantum 
gravity modifications of low-energy electrodynamics, and inverse quadratic dependence, mo- 
tivated by Faraday rotation generated by primordial magnetic fields. We constrain the 
parameters associated to each kind of dependence and use our results to give constraints on 
the models mentioned. We forecast the sensitivity that Planck data will be able to achieve 
in this respect. 

I. INTRODUCTION 

The improvement of experimental techniques for measuring the Cosmic Microwave Back- 
ground radiation (CMB) has allowed to measure with increasing accuracy the parameters of 
the standard cosmological model [1]. In particular, a boost to this field was given by the ability 
to not only measure the CMB temperature anisotropics, but also its polarization, whose under- 
lying physics is now very well understood, and which can give valuable information on crucial 
aspects of the universe history, such as the reionization epoch and the inflationary mechanism. 

The precision achievable with the last generation of experiments allows us to be ambitious 
and start considering CMB data as a tool to test also models of non-standard physics. In par- 
ticular, it is possible to constrain cosmological birefringence effects, and since the first WMAP 
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and BOOMERanG data were analyzed with this aim [2J, tens of papers appeared showing how 
very significant constraints can be already put on this kind of effect. 

Birefringence produces a rotation of the polarization direction of a linearly polarized electro- 
magnetic (EM) field during its propagation. For CMB radiation, this causes a peculiar mixing 
between the E and B modes of the polarization field. In particular, if we indicate with a the 
amount of rotation, the effect on CMB angular power spectra can be formalized as: 

Cf E = Cf E cos 2 (2a) + Cf B sin 2 (2a) 
Cf B = Cf E sin 2 (2a) + Cf B cos 2 (2a) 
C EB = \(c EE -C BB ) sin (4a) 
Cj E = Cj E cos (2a) 

Cf B = Cf E sin (2a) (1) 

where Cg refer to the spectra that one would expect without any birefringence effect, and Cg are 
the measured spectra if a polarization rotation has indeed happened. Note that in particular 
birefringence predicts non zero TB and EB cross-correlation, which would be instead vanishing 
according to standard parity-preserving physics. 

Motivations for testing cosmological birefringence arise in many contexts of theoretical 
physics. In general a mixing between E and B modes, which have opposite properties un- 
der parity transformations, signals the occurrence of some parity-violating phenomenon, but 
the details of the mechanism can be very different. 

Already standard physics predicts a mixing between E and B modes as described in eq. 0, 
if CMB photons happen to travel through a magnetic field. This effect is known as Faraday 
rotation, and can affect the whole CMB sky if caused by primordial magnetic fields OH]. Be- 
sides this, there are several non-standard physics models implying birefringent behavior of light, 
which are well motivated from a theoretical perspective. Some examples are the coupling of 
the electromagnetic (EM) field with a quintessence field [5-7J, coupling of the EM field with 
axions [8], Chern-Simons interactions in the EM Lagrangian [91-41 1 j and quantum-gravity moti- 
vated effective theories for electromagnetism, introducing non-renormalizable, parity- violating, 
corrections to the EM Lagrangian |12H14| . 

This plethora of different models calls for finding a way to discriminate between them if any 
birefringence signature in CMB polarization data appears. We will focus on the main feature 
that could make this possible, at least in principle, i.e. the energy dependence of the induced 
polarization rotation. 
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Presently available constraints, obtained by combining different CMB datasets, assume 
energy- independent birefringence (the most recent and complete analysis is (Ej). This can 
provide general information on the strength of the birefringence effect, but is of course very 
limiting for testing specific models predicting non-trivial energy dependencies. 

We will show that, although the quality of current CMB data is not yet good enough to 
discriminate among the energy dependencies investigated here, constraints on the relevant the- 
oretical models may be significantly improved by properly accounting for the observational 
frequency bands. We will also show that a great improvement is to be expected once Planck 
data will be analyzed in the way we suggest. 

The next section provides a review of current constraints on birefringence coming from 
the analysis of WMAP7, BOOMERanG, QUAD and BICEP data. We also give a review on 
the combined analyses proposed so far. In section III we review non-standard models which 
predict specific energy dependence of the birefringence effect. We then propose a method to 
combine datasets (different experiments and difference frequency channels) according to the 
energy dependencies we want to test and we provide our constraints. In section IV we forecast 
on the sensitivity reachable by the Planck mission, if an appropriate treatment of the information 
coming from different frequency channels is performed. We show that the sensitivity on energy- 
dependent effects will in general improve of about one order of magnitude with respect to what 
it is achievable with CMB present data. Section V provides a discussion of our results in light 
of the different models we have tested, and translates our results into constraints on the specific 
model parameters. We show that for some of the models our method allows to improve present 
constraints of about one order of magnitude. 

II. CURRENT CONSTRAINTS 

In this section we briefly report and discuss present available constraints on birefringence 
coming from CMB polarization observations. 

Since our main goal is to provide an energy dependence analysis, ideally we would only 
report results coming from the analysis of single energy channels. However, because not all the 
constraints available in the literature are of this kind, we will also consider those drawn from 
multifrequency experiments, where the different channels have been combined. The issue of 
accounting for the frequency dependence implicitly encoded in such results will be discussed in 
the next section. 
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The most recent constraints on a from CMB observation are shown in Table |IJ which also 
shows the frequencies observed, the multipole range explored and the uncertainty associated to 
the estimates. 

TABLE I. Presently available constraints on birefringence angle a from CMB. We report the experi- 
ment from which data were taken, the observational frequency channels, the harmonic multipole range 
explored, the estimated value of the birefringence angle, together with its statistical and systematic un- 
certainty (unless differently stated in the footnotes). The lines in boldface correspond to the data that 
we actually use for our analysis (see the text for further details). 



Experiment 


Energy channels (GHz) 


t 


a ± stat(±syst) (deg) 


Reference 


WMAP7 


33+41+61 


2-23 


-3.8 + 5.2 + 1.5 


DP 


WMAP7 


41+61+94 


24-800 


0.9 + 1.4+1.5 


m 


WMAP7 


33+41+61+94^] 


2-800 


-f. 1 + 1.4+1. 5 


m 


WMAP7 


33+41+61 


2-23 




m 


WMAP7 


33+41+61 


2-47 


-1.6 + 1.7 


on 


WMAP7 


33+41+61 


2-30 


A 9 + 1.9+10.2 
^• z -3. 1-7.5 


[13 


WMAP7 


33+41+61 


2-800 


i O+0.6+2.3 
x -°-0. 7-2.3 


on 


BOOM03 


145 


150-1000 


4.3 + 4.1^] 


m 


QUAD 


100 


200-2000 


1.89 + 2.24 + 0.5 


m 


QUAD 


150 


200-2000 


0.83 + 0.94 + 0.5 


m 


QUAD 


100+150 


200 - 2000 


0.64 + 0.5 + 0.5 




BICEP 


100+150 


21-335 


2.60 + 1.02 + 0.7 


m 





a The lower (higher) frequency does not contribute to the high (low) multipole constraint (cfr. the first two 
lines in the table) 

b For this result and the following one should add the same systematic as for other WMAP7 results 
c This error already takes into account BOOMERanG systematic error of —0.9° ± 0.7° 
d The systematic uncertainty is given in |21| 

For QUAD data we report both the combined-channel result and the constraints coming 
from the two single-channel analysis. The former results from an updated analysis, and gives 
in principle stronger constraints. Unfortunately, no update on the single frequency results has 
been reported by the QUAD team. Since we are mostly interested in analyzing potential energy 
dependence of the birefringence effect, we will in the following work with the QUAD single 
channel results, even if they are less recent. 

Constraints on a from WMAP7 data have been put either on the whole range of multipoles 
accessible, or by focusing on the low and high multipoles separately. Multipoles below i ~ 
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20 are expected to be mostly sensitive to photons that underwent reionization, at z ~ 10. 
Since reionization resets the polarization anisotropies produced at the last scattering surface, 
the observed birefringence effect would result as a superposition of contributions belonging to 
different scales and times. In order to make our dataset uniform and concentrate only on the 
energy dependence of birefringence, we will include in our study only results obtained for the 
multipole range 24 -j- 800. 

The WMAP team also reports on a combined constraint on the birefringence angle from 
WMAP7, QUAD and BICEP, where an inverse variance weighting scheme is adopted and 
systematic errors are taken into account by adding them (in quadrature) to the statistical ones 
[1]. They find a = —0.25° ± 0.58°, which is well compatible with no polarization rotation 
effect. In [15] the constraint from WMAP7, BOOMERanG and BICEP observations is a = 
—2.28° ± 1.02°, which results from a joint analysis of all angular power spectra, and takes into 
account also systematic uncertainties. This gives a 2.2 standard deviations indication for a 
birefringence effect. However, such indication is strongly reduced by the addition of QUAD 
data, providing a = —0.26° ± 0.54°, including systematic uncertainties. 

This is of course due to the fact that the estimate of the polarization rotation angle com- 
ing from QUAD data (combined energy channels) has some tension with results from other 
experiments (a ~ 2a tension under the assumption of energy-independent effect [IS]), which 
comes mostly from the 150 GHz channel (see Table [I]). A brief discussion about such tension is 
provided in [15] and [llj^| We will come back on the topic in the next section. 

All the constraints discussed above, obtained combining data from different experiments, 
assume no energy dependence of the birefringence rotation angle. This is of course limiting if 
one wants to put reliable constraints on an effect which might very well have some non trivial 
dependence on the radiation energy. In the next Section we will focus on how to exploit current 
estimates of a to investigate its energy dependence. 

III. STUDY OF POSSIBLE ENERGY DEPENDENCE OF BIREFRINGENCE 

EFFECTS 

In this section we test the hypothesis that birefringence of CMB photons might have some 

specific energy dependence. As we mentioned in the introduction, this can be a powerful method 

1 Note that in |15j also constraints from the single experiments are re-derived, both including and not including 
systematic errors. These are in very good agreement with the ones reported in our Table [I] once statistical 
and systematic errors are simply summed in quadrature. 
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to discriminate between different theoretical models predicting a rotation of polarization during 
photons propagation. Here we will limit ourselves to qualitatively describe the dependence of 
the polarization rotation angle a as a function of photon energy E. In the final discussion 
section we will make actual connection between these qualitative behaviors and specific models, 
relating the parameters used in this and the following section with the parameters characterizing 
the models considered. 

The laws of energy dependence that we will test are of four kinds: 

• square law dependence (ot(E) = qE 2 ). This is predicted by quantum gravity-inspired 
models, such as the Lorentz violating effective field theory for Electrodynamics proposed 
in [12] and was already considered in a cosmological framework in [13^. E] . 

• linear dependence (a(E) = IE). This kind of dependence is motivated by quintessence 
coupling models (eventually string-inspired) [6], and more in general by Lorentz violating 
electrodynamics, such as renormalizable 'Weyl' interactions [131 122] . 

• no dependence (a(E) = ao). This is the most studied case in the CMB literature, 
and basically the only one considered when different datasets are combined (with a few 
exceptions, such as a work by one of the authors [13]) Energy independent birefringence 
is expected in the presence of Chern-Simons interactions of the EM field [9| 122]. 

• inverse square law dependence (ot(E) = hE~ 2 ). This is the kind of dependence expected 
from Faraday rotation due to the photons traveling through a magnetic field. In partic- 
ular a uniform rotation across the sky can be traced back to the presence of primordial 
homogeneous magnetic fields. 

A. Combining different datasets 

Our approach will be to fit for the four kinds of energy dependence described above using the 
available constrains on birefringence angle coming from individual CMB experiments, reported 
in boldface in Table HI 

In combining the constraints coming from different experiments in such a way that allows 
for a study of the energy dependence we have to face two issues. The first one is that the each 
energy channel to which a CMB experiment is sensitive will have a certain width. Instrumental 
response per frequency is usually very well modeled by a constant function within a given 
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frequency range, its bandwidth, which is of order 20 -j- 30% of the nominal energy. In Tab{IT|we 
report the frequency channels used by the various experiments to constrain birefringence and 
the channels bandwidth. The second issue is the fact that for some experiments (WMAP7 and 
BICEP, see table [I]) the constraints are derived from a combination of several channels, under 
the assumption that the amount of polarization rotation does not depend on energy. 

TABLE II. Energy channels and corresponding bandwidth of the experiments considered in the present 
work. 



Experiment 


Energy (GHz) 


Bandwidth (GHz) 


WMAP7 


41 


8 




61 


14 




94 


20 


BOOM03 


145 


45 


QUAD 


95 


26 


QUAD 


150 


41 


BICEP 


96 


22 




150 


39 



Let us formalize the issue of properly accounting for multifrequency and band widths. If the 
true functional dependence of a on energy is a = a(E), then we can interpret the measured 
value af from each experiment as the result we would have obtained if we had only considered 
data at some effective energy Ef, i.e. at = ct(Ef). Of course the actual value of Ef depends on 
the specific functional dependence of a on E. The relation between ctf and the function a(E) 
can be written as: 

a f = ^J a(E)f(E)dE, N = J f(E)dE (2) 

where f{E) is a weighting function. We assume a constant response for each frequency channel i 
within a given energy interval [E_ , E^] (its bandwidth) and we will adopt a weighting scheme 
of the different channels according to the corresponding pixel nominal sensitivity [^] So: 

, I la} if Ee [E^,e¥ ] ] 
f(E) = { ' 1 (3) 
elsewhere 
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This is coherent with the assumption, made in [T], that the uncertainty on a is related to the signal to noise 



ratio of TE and EE in the following way: Err[a TB ] ~ 1 TE , Err[a EB] ~ 



2(S/JV) J * > ' ' L~ ' — 2(S/N) 1 
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where Oi is the instrumental noise per pixel. From Eqs. Q and ([3]), the effective energy Ej 
associated to an experiment to can be computed for any kind of assumed energy dependence. 

For clearness, we provide in Appendix [A] the explicit example of a multi-channel experiment 
(as WMAP and BICEP for the purposes of the present work), for the different kinds of energy 
dependence under investigation. 

By means of the procedure sketched above, for each functional energy dependence we want 
to test we translate the information of Table [TT] into an effective energy Ej to be associated to 



the value of a measured by each experiment. Results are reported in Tab. Ill, together with 
the corresponding constraint on a and the associated error. The latter includes statistical and 
systematic uncertainties, taken as independent and so summed in quadrature. 

To test a specific energy dependence, we simply perform a best-fit analysis of the data in 



Table III to constrain the parameter of the models, o>q, I, q, h. 



TABLE III. For each experiment considered, we report the effective energy (for the different kinds of func- 
tional energy dependence under investigation, except the energy independent case), the corresponding 
measure of a and the associated error (statistical and systematic combined). 



Experiment 


Energy (GHz) [lin/quad/inv] 


a (deg) 


a(a) (deg) 


WMAP 7 


53.1 / 55.4 / 47.7 


-0.9 


2.0 


B03 


145.0 / 147.3 / 137.8 


-4.3 


4.1 


BICEP 


135.9 / 139.5 / 124.0 


-2.6 


1.2 


QUAD1 


95.0 / 96.2 / 91.3 


-1.9 


2.3 


QUAD2 


150.0 / 151.9 / 144.3 


0.8 


1.1 



We will treat each experiment as independent. In fact, each instrument is characterized 
by a different sky coverage and instrumental noise properties. Confirmation of the fact that 
the dataset can be safely assumed as independent comes from the agreement between the two 
combined analysis of |15| and pQ. The former performs a joint analysis of all the angular power 
spectra, while the latter only combines the final results derived from each experiment. 

Note that for QUAD, that provided a separate analysis of the two frequency channels, we 
consider the two data independently, to improve tests of energy dependence. Moreover this 
allows us to study the influence of the 150GHz channel data on the final fits. As we have 
mentioned in the previous section, this channel provides a constraint on a that has a ~ 2a 
tension with the combination of all other data, both from other experiments and from the 
100GHz channel of the same experiment. A detailed analysis of the power spectra arouses 
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suspicions about the possible presence of unaccounted systematic errors [IT]. So we think 
it is worth to treat this issue in a careful way, to avoid hasty conclusions about the fate of 
birefringence models. 



B. Results 



Our best-fit parameters are reported in Table IV, for two cases: including (left) or not 
(right) the QUAD 150GHz channel in the analysis. We have fitted for energy-independent 
behavior (a(E) = ckq), linear energy dependence (a(E) = IE), quadratic energy dependence 
(a(E) = qE 2 ) and inverse quadratic dependence (a(E) = hE 2 ). Errors on fit parameters have 
been computed as shown in Appendix |B} 





with QUaD 150 GHz 


without QUaD 150 GHz 


parameter 


fit value ± error 


reduced x 2 


fit value ± error 


reduced \ 2 


a [deg] 


-0.93 ±0.70 


1.4 


-2.22 ±0.92 


0.17 


/ [deg -GHz- 1 ] 


-6.1 • 10~ 3 ± 5.3 • 10~ 3 


1.5 


-1.99 -10- 2 ± 0.78 -10- 2 


0.03 


q [deg • GHz" 2 ] 


-3.6- 10~ 5 ±3.6- 10~ 5 


1.6 


-1.44- 10~ 4 ± 0.57 • 10~ 4 


0.06 


h [deg • GHz 2 ] 


-4.2- 10 3 ±4.2- 10 3 


1.6 


-5.0- 10 3 ±4.3- 10 3 


1.32 



TABLE IV. Fit parameters, together with their uncertainties and corresponding reduced \ 2 ■ On the left 
results obtained by including QUaD 150 GHz channel in the analysis and on the right results obtained 
not including that channel. 



Our best-fit results are also shown in Figs. [T]and[2j together with data points and associated 



errors (taken from Tab III). Two curves are shown in each panel: in (dashed) blue the best-fit 
curve obtained by including the QUAD 150 GHz channel (blue point) in the dataset; in (solid) 
black, the best-fit curve obtained after rejecting the QUAD 150 GHz channel. The black thin 
line is for the reference value a = 0. 

We see that in general, as expected, excluding the QUAD 150 GHz channel gives constraints 
on the parameters which are more significantly different from zero. In this case, while the linear 
and quadratic energy dependence parameters are both at a bit more than 2.5<r from zero, the 
energy-independent case parameter is at 2.4<r and the inverse energy dependence parameter is 
at only l.lcr from zero. We can not say anything significant about the preferred model, since 
the x 2 values are all very small, except for the inverse quadratic case, which is also the one 
whose parameter is estimated to be closer to zero. Nevertheless, looking at the plots we see that 




20 40 60 80 100 120 140 160 20 40 60 80 100 120 140 160 

Energy [GHz] Energy [GHz] 



FIG. 1. Left panel: no energy dependence. Right panel: linear energy dependence. In both figures the 
thin black line marks the zero, the (dashed) blue line is the best fit curve obtained including the QUAD 
150 GHz channel (blue point), while the thick black line is the one not including the QUAD 150 GHz 
channel. 




20 40 60 80 100 120 140 160 20 40 60 80 100 120 140 160 

Energy [GHz] Energy [GHz] 



FIG. 2. Left panel: quadratic energy dependence. Right panel: inverse quadratic dependence. For 
details see caption of Fig. [I] 



indeed, if one considers only the black points, they appear to follow very well both the linear 
and quadratic dependence black curves (better than the other kinds of trend), but their error 
bars are too big to make any serious claim. We would only take this observation as an incentive 
to look at energy dependence of future experiments data, which are expected to provide a 
single-channel sensitivity on a about one order of magnitude better (see next section). In that 
case, it will hopefully be possible to make claims about the trend of the polarization rotation 
with respect to energy. 

Including the QUAD 150GHz channel in the analysis gives on the other hand reduced x 2 
values around 1.5, and fit parameters closer to zero: slightly more than la from zero for the 

10 



constant and linear case, la for the quadratic and inverse quadratic case. 

We can also estimate the probability that available data are compatible with no effect at 
all, regardless of a possible energy dependence. To this aim we sum the residuals of each data 
point from the zero line (xfj) an d provide the corresponding probability. Including all QUAD 
points we get Xo = 7-27 (probability 20%), while excluding the 150GHz channel gives Xo = 6.67 
(probability 15%). In both cases data disfavor the hypothesis that what we see is only due to 
fluctuations around zero, but we can not exclude it yet. 

In the conclusions we will translate our findings into constraints on parameters of the models 
discussed at the beginning of this section. Here we just want to stress how crucial it seems to 
understand the reliability of the QUAD 150GHz channel one one hand, and on the other to 
improve the sensitivity to birefringence for the single energy channels of future experiments. In 
fact it seems that the possibility to significantly discriminate between different models, giving 
stringent constraints on their parameters, is very close. 

IV. PLANCK FORECAST 

In this section we forecast the absolute error on fit parameters for the Planck mission |23j . 
We focus on the central frequency bands explored by the Planck satellite, 70GHz, 100GHz and 
143GHz, whose frequency bandwidths are, respectively, 14GHz, 33GHz and 48GHz. We refer to 
|14j for the sensitivity to the birefringence angle for each channel. We follow the same procedure 
adopted for the other CMB experiments to account for the observational bandwidth, keeping 
the three energy channels separated in order to have a wide coverage in the frequency domain. 
The error bars associated to the best fit parameters are estimated exactly in the same way as 
we have done above (see Appendix [B|) , and compared to the previous results. 

TABLE V. For each Planck channel considered, we report the effective energy (for the different kinds 
of energy dependence under investigation, except of course for the energy independent case), and the 
corresponding sensitivity on a 



Planck channel (GHz) 


Effective Energy (GHz) [lin/quad/inv] 


cr(a) (deg) 


70 


70.0 / 70.5 / 68.6 


0.64 


100 


100.0 / 101.8 / 94.4 


0.14 


143 


143.0 / 145.6 / 135.1 


0.073 
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TABLE VI. Planck forecasts on fit parameter errors. Current sensitivities are based on the full dataset 
of Tab |III| Planck sensitivities are based on Planck instrumental properties as reported in TabfV] 



parameter 


Current sensitivity 


Planck sensitivity 


<j(a ) [deg] 


7.0 • 10- 1 


6.4- 10~ 2 


a(l) [deg ■ GHz^ 1 ] 


5.3 • 10~ 3 


4.8 ■ 10~ 4 


a(q) [deg • GHz~ 2 ] 


3.6 • 10~ 5 


3.3 • 10~ 6 


(r(h) [deg • GHz 2 ] 


4.2 • 10 3 


8.7- 10 2 



In Tab. (VI) we report both the current sensitivity on fit parameters and the forecasted 
sensitivity for Planck. 

We see that using Planck data by themselves will allow to improve all constraints of about 
one order of magnitude. The availability of several distinct energy channels will also allow to 
use energy dependence to discriminate between the different energy dependence possibilities if 
any kind of effect is found. 



V. DISCUSSION AND CONCLUSIONS 



In this work we have fitted available data on birefringence coming from CMB polarization 
measurements against several possible relations between the birefringence effect and photons 
energy. This is the first time that the issue of energy dependence is treated in such a systematic 
way. 

In principle this kind of analysis would allow to disentangle effects of different origin that 
result in the same kind of rotation of the polarization direction (like Chern-Simons interactions 
of the EM field, propagation of CMB photons through primordial magnetic fields, Planck scale 
Lorentz violations just to mention a few). Even if polarization data are not good enough to prefer 
a specific trend among the competing ones investigated in this work, we are able nonetheless 
to improve constraints on the underlying theoretical models up to one order of magnitude with 
respect to previous results. 

For example, we have studied the possibility of a quadratic energy dependence of the bire- 
fringence effect. This might be traced back to Quantum Gravity Planck-scale effects [12} I14j. 
whose relics at low energies can be modeled as a coupling of the EM field with a fixed time- 
like vector. If, following |12j and [Hj, we write the coupling constant between the EM field 
and the vector through a dimensionless parameter £ and the Planck mass scale Mp, then the 
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polarization rotation angle is related to the distance covered by photons and their energy by: 

a(E) = -^-E 2 Al (4) 
y ' M P y J 

Taking into account photon energy redshift from the last scattering surface (zls) this becomes 

a(E ) = jj-^o Jo 2iS ( 1 + z)H(z)~ 1 dz, where E is the photon energy today. So our constraint 

q = (-3.6±3.6)-10" 5 deg-GHz- 2 can be translated^ into a constraint on the Quantum Gravity 

parameter £ = —0.22 ± 0.22. If we disregard the QUAD 150GHz channel, the constraint 

q = (-1.44 ± 0.57) • 10 _4 deg • GHz~ 2 translates into £ = -0.87 ± 0.34. To compare with 

previous constraints from CMB [14J, £' = —0.110 ± 0.075, we have to take into account the 

different convention used in treating the Planck mass scale. In [14J it was used the reduced 

Planck mass Mp^ = Mp/^/8tt. We see in this way that our constraint improves the previous 

one of about a factor three. 

Concerning the linear energy dependence case, our results can be translated into constraints 

on the 'Weyl' interaction described in [22] (see also [13]). We are not giving here the details of 

the model, which is well discussed in the two references given above. We just report that, in 

this case, the polarization rotation angle is related to the distance Al traveled by photons and 

the dimensionless scalar ^o> which sets the amplitude of the interaction, by 

\a(u)\ =4u\V \A£. (5) 

Taking into account the Universe expansion and the energy redshift, our parameter for the linear 
dependence is related to through |/| = 87r| ^ o I Jq LS H(z)~ 1 dz. Assuming gaussian errors for 
our parameter I, so that the limit on its absolute value is \l\ < 8.7 • 10 _3 deg/GHz at 68% c.l. 

< 1.5-10- 2 deg/GHz at 95% c.l.), we get |* | < 5.8- 10~ 33 ft, at 68% c.l. (|tf | < 1.0-10~ 32 /io 
at 95% c.l. ), where ho is the reduced Hubble constant^} Such a constraint improves the previous 
one |13] of about one order of magnitude. Disregarding the QUAD 150GHz channel, our 68% 
c.l. upper limit on I is |/| < 2.4 • 10" 2 deg/GHz (|Z| < 3.3 • 10~ 2 deg/GHz at 95% c.l.), so 
|^o| < 1.6 • 10~ 32 ho at 68% c.l. (|* | < 2.2 • 10~ 32 h at 95% c.l.) 

The constraint on the inverse quadratic energy dependence of the birefringence effect can be 
translated into a constraint on the amplitude of uniform primordial magnetic fields. Using the 
result of [3] and [24] . we can relate the rms value of the rotation angle to the intensity of the 
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Here and in the following we assume a cosmological model with Ho = ho ■ lOOKm/s/Mpc = 71 Km/s/Mpc, 



ft m = 0.27 and ft A = 0.73. 

4 Note that this kind of 'Weyl' interaction would also produce a leakage from linear to circular polarization, 



which could be detectable by looking for a circular component of the CMB. 
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primordial magnetic field through: 

<« 2 > 1/ ^ 1 - 3 °( 1 0^)(30ik)" W 

So our parameter /i is related to the magnetic field intensity in the following way: < h? > x / 2 = 
(icF^Gauss) (30GHz) 2 . From our constraint on h, assuming gaussian errors we derive the 
upper limit < h 2 > 1 t 2 < 6.3- 10 3 at 68% c.L, which translates into the constraint on the magnetic 
field intensity Bq < 5.4- 10 -9 Gauss and Bq < 6.1 • 10~ 9 Gauss at 68% c.L, with and without the 
QUAD 150GHz channel respectively. Of course this result is valid only under the assumption 
of a uniform primordial magnetic field. In most general scenarios the magnetic field will have 
some non trivial spectral distribution, which would result in a direction-dependent birefringence 
(see e.g. |25j). 

As last case we consider energy independent birefringence, which is the most broadly studied 
in the literature. We can translate our constraint ao = (—0.93 db 0.70)°, obtained including the 
QUAD 150GHz channel in the analysis, into an upper limit \a\ < 1.26° at 68% c.L and 
\a\ < 2.11° at 95% c.L This can be used to constrain the Chern-Simons interaction of the 
EM field described in [9]. The relevant parameter of the model is \Lq — Lcos<j>\, which sets 
the intensity of the interaction. The birefringence angle is related to this parameter and the 
distance covered by photons (A£) by 

\a\ = ^\L -Lcos<P\A£^ (7) 

Taking A£ = ^ LS {l+z)- 1 H(z)- 1 dz = 2.15-KT 33 /i eV we find \Lq-L cos <f>\ < 9.4-10~ 44 /io GeV 
at 68% c.L and \Lq — L cos <p\ < 1.6 • 10 _43 /io GeV at 95% c.L This constraint is slightly stronger 
than the one found in [13], since we are taking into account a larger dataset. When disregarding 
the 150GHz QUAD channel, our constraint ao = (— 2.22 d= 0.92)° translates into the upper limit 
\a \ < 2.66° at 68% c.L, and \a \ < 3.77° at 95% c.L In this case |L o -Lcos0| < 2.0-10" 43 /i GeV 
at 68% c.L and \L - Lcos(p\ < 2.8 • 10~ 43 /i GeV at 95% c.L Note that limits on the constant 
energy dependence have to be taken with a grain of salt, since the same kind of polarization 
rotation might also be mimicked by systematic errors in polarimeters calibration [18J. 



<j> is the angle between a preferred direction pointed out by the vector and the photon propagation direction. 
Since we use data averaged out over the whole sky, the limits we put have to be interpreted as limits on the 
average value of \Lq — Lcos<j)\ over all directions. 
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Appendix A: Multi-channel experiments 

In this appendix we provide explicit calculations of the effective energy for a multi-channel 
experiment. The four energy dependence of a under investigation are considered. 

• No dependence (a(E) = ao). This case is trivial, since aj = ao and no effective energy 
is introduced. 



• Linear dependence (ct(E) = IE). 

I lEf(E)dE 
° f= Jf(E)dE = lE f^ E f 

• Quadratic dependence (a(E) = qE 2 ). 

_ fqE 2 f(E)dE _ 2 
af ~ Jf(E)dE =« E f^ E f 

• Inverse square law (a(E) = hE~ 2 ). 
JhE- 2 f(E)dE 



E< h (4 4) - E - 



i \ 



(Al) 



3 Z^(e?-e 



(0 



1/2 



(A2) 



a f 



Jf(E)dE = kE f * E ' 



Ei ^ V 1 - i E +)~ 1 

i * 



£ W _ e ® 



-1/2 



(A3) 



For WMAP, the pixel noise level of the channels we consider is [22]; o"i = 2.672 mK for the 
Q-band (41 GHz), a 2 = 3.371 mK for the V-band (61 GHz) and a 3 = 6.877 mK for the W-band 
(94 GHz), while for BICEP o\ = 0.81 fiK for the 96 GHz channel and a 2 = 0.64 /M for the 150 
GHz channel [27]. 

Appendix B: Fit parameter errors 

In this appendix we provide how to compute the error associated to the fit parameter for the 
four energy dependences considered in the present work. Here Xi is the i-th value (data point) 
we fit and Oi is the associated error. 



15 



No dependence (a(E) = ao). 

'(«)-(x:^) 

Linear dependence (ot(E) = IE). 

*(0 = f £ f 

Quadratic dependence (a(E) = qE 2 ). 



Inverse square law (a(E) = HE 



-1/2 



1/2 



07 



*(<?)=£ 



-1/2 



-1/2 



(Bl) 



(B2) 



(B3) 



(B4) 
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